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[57] ABSTRACT 

By using two SAR antennas spaced a known distance, 
B, and oriented at substantially the same look angle to 
illuminate the same target area, pixel data from the two 
antennas may be compared in phase to determine a 
difference A<|) from which a slant angle 6 is determined 
for each pixel point from an equation A</>=(27rB/\)sin(- 
0 — a), where X is the radar wavelength and a is the roll 
angle of the aircraft. The height, h, of each pixel point 
from the aircraft is determined from the equation h = R 
cos 0, and from the known altitude, a, of the aircraft 
above sea level, the altitude (elevation), a', of each point 
is determined from the difference a — h. This elevation 
data may be displayed with the SAR image by, for 
example, quantizing the elevation at increments of 100 
feet starting at sea level, and color coding pixels of the 
same quantized elevation. The distance, d, of each pixel 
from the ground track of the aircraft used for the dis- 
play may be determined more accurately from the equa- 
tion d=R sin 6. 

10 Claims, 4 Drawing Figures 
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METHOD AND APPARATUS FOR CONTOUR 
MAPPING USING SYNTHETIC APERTURE 
RADAR 
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ORIGIN OF INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Section 305 of the National 
Aeronautics and Space Act of 1958, Public Law 85-586 10 
(72 Stat. 435; 42 USC 2457). 

BACKGROUND OF THE INVENTION 

This invention relates to synthetic aperture radar 
(SAR), and more particularly to a method and appara- 15 
tus for determining the elevation of terrain features in 
SAR image data, and if desired for determining more 
accurately the distance of each pixel from the ground 
track of the SAR platform. 

The development of SAR has made it possible to 20 
quickly map vast areas for geological as well as for 
defense purposes. Since topographic maps are much 
more useful, it would be important to use SAR for 
contour mapping. However, in mapping rough terrain, 
there is a problem in that terrain features of high alti- 25 
tude tend to be foreshortened (i.e., appear to be closer 
to the ground track of the SAR platform than they 
actually are), and the greater the slope of the terrain 
feature, the more pronounced the foreshortening effect. 
This is due to the fact that the range of each image pixel 30 
from the flight path is essentially determined by the 
slant range of the SAR platform to the reflecting points 
on the terrain. 

It would be desirable to determine the altitude of 
each point on the slope, and to determine the true hori- 35 
zontal distance of each point from the ground flight 
path, in order to present an image that accurately por- 
trays the contours in rough terrain without foreshorten- 
ing. 

40 

SUMMARY OF THE INVENTION 

In accordance with the present invention, two SAR 
antennas are spaced apart a distance, B, on a platform 
which carries the SAR over terrain at a known altitude, 
a, with respect to a fixed reference, such as sea level. 45 
The return signals from each point to the antennas are 
separately processed to the pixel level, and the phase 
differences in the signals returned through the two 
antennae from substantially the same points are de- 
tected. The phase difference, A<J>, detected for each 50 
point is related to the spacing B between the antennas 
by the following equations: 

. , 2ttB . . 

A© = " sin ( 6 — a) 55 

where 0 is the angle between the slant range line and a 
vertical line, a is the roll angle of the SAR platform and 
X is the radar wavelength. The height of the SAR plat- 
form above each point is related to the slant range by 60 
the equation 

h—R cos 0 

Thus, from the measured phase difference, A</>, the slant 65 
angle 0 is determined for each pixel from the first equa- 
tion. Then the height, h, of the platform is determined 
from the second equation. Since the altitude, a, of the 
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SAR platform is known from an independent instru- 
ment measurement, the altitude, a', of each pixel (point) 
is determined directly from the equation 

a’ = a—h. 

For more accurate position display of the pixels, the 
distance, d, of each point is determined from the equa- 
tion 

d—R sin 6 

Each more accurate pixel may then be properly dis- 
played in distance from the ground track of the flight 
path with its altitude indicated. Contour lines may be 
displayed by connecting points of equal altitude a' at 
regular quantized increments of altitude (elevation), or 
by simply color coding pixels at regular quantized in- 
crements of elevation. 

BRIEF DESCRIPTION OF DRAWINGS 

FIG. 1 is a diagram useful in understanding the nature 
of the problem which the present invention presents. 

FIG. 2 is a block diagram of an exemplary implemen- 
tation of a SAR for the present invention. 

FIG. 3 is a block diagram of an off-line system for 
processing image data in accordance with the present 
invention. 

FIG. 4 is a geometric diagram useful in understanding 
the technique of the present invention, which is a 
method for determining true range (or distance from the 
ground track of the SAR flight path) and altitude of 
each pixel. 

DESCRIPTION OF PREFERRED 
EMBODIMENTS 

Referring now to the drawings, FIG. 1 shows the 
geometry involved in the foreshortening effect of a 
SAR system carried by an aircraft 10 over rough terrain 
12 illustrated as a slope rising from sea level, or some 
other known reference level. Doppler processing local- 
izes a given pixel to the plane of the figure, and range 
processing localizes each pixel to an arc, such as pixels 
at points 13 and 14 part of the way up the slope which 
are localized to arc 15 and 16. 

It is clear that the point 14 has a slant range, R, 
greater than the slant range R' of point 13. Since the 
slant range R for the point 14 is only slightly greater 
than at the point 13, it will be displayed at a distance 
from the ground track of the aircraft less than the actual 
distance d from the ground track, but greater than the 
distance of point 13. This foreshortening effect may be 
corrected as a consequence of the present invention by 
which the height, h, of the SAR above the point 14 is 
determined, and from its height and slant range, both 
the altitude a' and the distance d of the point 14 from the 
ground track are determined. The pixel from point 14 
can be displayed with its true altitude a' at the proper 
distance from the ground track. 

The object of the invention is to determine the alti- 
tude a' of each point. As just noted, it is necessary to 
first determine the height, h, of the aircraft above the 
point, and from the known altitude of the aircraft deter- 
mine altitude a' as the difference between a and h. This 
is a geometric problem that can be solved with a pro- 
grammed SAR data processor once the slant angle 0 of 
the line of sight from the aircraft to the point is deter- 
mined. Note that this angle is distinct from the “look 
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angle,” which is the angle that the bore axis of the side- 
looking SAR is depressed below a plane passing 
through the wings of the aircraft. The look angle is 
fixed, but the slant angle of signal return from a given 
point at a particular range is a function of the height of 5 
the aircraft with respect to that given point. Conse- 
quently, the problem is to determine the slant angle 6 
for each point from which the height, h, and the dis- 
tance, d, of that point can be determined for displaying 
each pixel in its proper position with respect to the 10 
ground track of the aircraft, and to plot or otherwise 
indicate its altitude. 

To determine the slant angle 6 , a planar array antenna 
21 is connected to a standard SAR transmitter and re- 
ceiver, and a second planar array antenna 22 spaced a 15 
distance B from the first is connected to the same trans- 
mitter and either to the same receiver through a multi- 
plexer, or preferably to a separate receiver. 

The principle involved is to illuminate the target area 
with radar pulses and receive echoes through the two 20 
antennae fixed on the aircraft at substantially the same 
look angle. The arrival time at the antennas of the ech- 
oes from a reflecting point would be slightly different 
because of the spacing, B, between the antennas such 
that the range to each reflecting point is slightly differ- 25 
ent. By maintaining that phase information through the 
whole receiving and recording process, the elevation a 
of each pixel in the radar image can be computed. 

In an exemplary implementation shown in FIG. 2, 
antennae 21 and 22 are multiplexed by an antenna selec- 30 
tor switch 23 and coupling circulators 23a and 23 b to 
sample the echo signal first from one antenna and then 
the other every other transmitted pulse. The pulse repe- 
tition frequency ranges from 800 to 1000 pps depending 
on aircraft ground speed which is sufficiently high for 35 
echo returns from successive transmitted pulses to be 
from substantially the same points. A timing unit 24 
receives a ground speed signal, V, from aircraft instru- 
mentation to control the frequency of a timing oscillator 
which controls the generation of PRF trigger pulses to 40 
enable a chirp generator 25 to transmit a burst of cycles 
from the stable oscillator 26 the phase of which is con- 
trolled by the timing unit 24. A modulator 27 is con- 
trolled by PRF trigger pulses to provide in a traveling 
wave amplifier 28 linear frequency modulation of the 45 
signal from the oscillator 26 gated by the chirp genera- 
tor 25. This FM modulation (chirp) allows for a 
matched filter in the data processor (FIG. 3) to produce 
a compressed pulse from the chirp for better range 
resolution in a customary manner. The PRF rate is 50 
slaved to the aircraft ground speed so that the target is 
illuminated with a substantially constant number of 
pulses per kilometer. 

The received echoes are passed through a receiver 
amplifier 29, and a mixer 30 which receives an appropri- 55 
ate frequency from the phase controlled oscillator 26 
via a frequency multiplier 31. The received echoes are 
downconverted to video in the mixer, amplified in a 
video amplifier 32 and digitized in vector form (real and 
imaginary components) in an analog-to-digital con- 60 
verter and formatter 33. The analog echoes thus con- 
verted to digital vectors are then stored in a high den- 
sity magnetic tape recorder 34 as two records for later 
processing. 

At the end of the data run, the tape is removed and 65 
interfaced with a playback unit 35 having a digital cor- 
relator 36, as shown in FIG. 3. The correlator performs 
a two dimensional compression which then allows ex- 
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traction from the vectors in digital form not only of 
phase, but also of the intensity of each picture element 
in the two records. These intensity values are then as- 
sembled by an image processor 37, into an image for 
each record. By consuming the images in a display unit 
38, it is possible to determine whether or not the two 
records are registered. If so, the two images reinforce 
each other for a clear and brighter image, if not, the 
result is loss of the image, and one of the records must 
be adjusted relative to the other until a clear and 
brighter image appears, i.e., until there is registration 
between the records on a pixel by pixel basis. 

Once registration is achieved, corresponding vectors 
of the pixels are compared by the image processor to 
determine the phase difference, A<f), between the anten- 
nae. The differential phase thus determined for each 
pixel is indicative of the height, h, of the aircraft above 
the pixel. Then the image processor determines for each 
pixel the altitude (elevation), a', and distance, d, of the 
reflecting point and encodes the pixel in accordance 
with both the distance, d, and the altitude, a', for display 
through the display unit 38 with every pixel at its 
proper distance and with altitude indicated, as by color 
coding the intensity display of points after quantizing 
their altitude at increments of 100 feet. 

From the foregoing, it is evident that the systems 
shown in FIG. 2 is composed mainly of standard radar 
circuit components. The features that are different from 
prior art imaging radar are the use of two spaced planar 
array antennas mounted so that they can receive echoes 
simultaneously, or substantially so, from the same target 
area. The antennas are mounted at a known fixed sepa- 
ration, B, on the aircraft for altitude (elevation) compu- 
tation by the image processor in the following manner. 
From the phase difference, A<f>, determined by compar- 
ing the vectors of the two echoes received through the 
two antennae, the slant angle 0 is computed from the 
following equation: 


A<f> = — ^ — sin (9 — a) 


( 1 ) 


where \ is the radar wavelength and a is the roll angle 
of the aircraft, as shown in FIG. 4. From that, the 
height of the aircraft above the reflecting point is com- 
puted from the equation: 

h = R cos e (2) 

And from the known altitude, a, of the aircraft, the 

altitude (elevation), a', of the point is computed from 
the equation: 

a‘=a—h (3) 

The distance d of the point from the ground track of the 
aircraft may be computed also from the following equa- 
tion: 

d—R sin 9 (4) 

This information (a' and d) is stored with the intensity 

value for each pixel. The value d may be used in plot- 
ting or otherwise displaying each element of the image 
through the display unit 38. 

In displaying the SAR image data, elevation may be 
quantized at intervals of 100 feet and displayed in differ- 
ent colors in order to provide one type of contour map, 
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as noted hereinbefore. To produce an image with more 
conventional contour lines, it would be possible to dis- 
play all elements at their proper distance and then su- 
perimpose lines connecting adjacent elements of the 
same quantized elements at intervals of 100 feet starting 5 
with the lowest elevation. Another way to display ele- 
vation information for quantized elements would be to 
print the quantized elevation at the center of every 
predetermined number of adjacent elements of the same 
quantized elevation. Still other ways of displaying ele- 10 
vation data will occur to those skilled in the art once the 
values d and a' are computed in accordance with the 
present invention for each element (pixel) by the image 
processor for display. In that manner, a contour map of 
vast areas may be produced in a relatively easy and 15 
quick manner with an accuracy which largely depends 
on the spacing, B, between the antennae. The greater 
the spacing, the greater the accuracy. 

FIG. 4 illustrates the problem in geometry which the 
image processor 37 computes for each pair of received 20 
echoes after having their phase difference Acf> and slant 
range R computed by the correlator 36 . The roll angle, 
a, taken into account in Equation (1 ) is exaggerated for 
emphasis in FIG. 4 . In practice, that angle will be very 
small, but since it is a parameter continuously moni- 25 
tored by the aircraft flight instruments, it is recorded 
with the SAR data, along with aircraft altitude so that 
roll angle may be taken into account for greater accu- 
racy. 

From Equation (1), the slant angle 0 is computed. 30 
Then using the slant range, R, the height, h, of the 
aircraft from the point P is computed from Equation (2). 
The altitude (elevation) a' of the point P is then com- 
puted as the difference a— h, as indicated by Equation 
(3). This provides elevation data for imaging, but to 35 
improve the accuracy of the contour image that can be 
produced from that image data, the correct distance, d, 
of each element is computed from Equation (4) for use 
in generating the image display. All of these equations 
are relatively easy to compute in an image processor 40 
comprised of a programmed digital computer. A simple 
technique would be to use a table for each equation. 

In an actual reduction to practice, the phased array 
antennae were placed symmetrically about the fuselage 
on the underside of the wing of a Convair CV-990 air- 45 
craft. The antennae were spaced 11.17 meters apart 
with the boresight of the antenna 21 on the right side 
(looking at the aircraft from the front at 24.5° from the 
vertical and the boresight of the antenna 22 on the left 
side 24° from vertical. The 0.5° boresight offset between 50 
the antennae minimizes the possibility of sidelobe inter- 
ference signals between the two antennae without de- 
grading the accuracy of the system for aircraft altitudes 
of more than about 2,000 and 3,000 feet above the ter- 
rain. 55 

Approximately twenty seconds of imaging data ob- 
tained over the Monterey Peninsula were processed 
into two images in complex format, one image for the 
data from each antenna. Then the pixel altitude a' was 
computed for each element and displayed in a color 60 
code in a manner described hereinbefore, starting with 
the lowest quantized altitude. (Image brightness for 
each pixel correspond to received echo power.) Blue 
was used for sea level, and increasing altitudes were 
shown in green, yellow, orange, magenta, purple and 65 
violet. The full band was thus equal to about 500 feet of 
elevation difference. The colors were then repeated. 
The color contours generally followed the local con- 
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tours with excellent agreement, even without comput- 
ing the distance d to correct the position of each pixel. 
Altitude variations of 50 feet are thus easily discerned. 
Sensitivity was greatest in the near range. 

Although preferred embodiments of the invention 
have been described and illustrated herein, it is recog- 
nized that modifications, improvements and equivalents 
may readily occur to those skilled in this art. For exam- 
ple, since several echoes are returned in succession from 
essentially the same point, pixels from the two records 
may be combined for the same point to achieve the best 
estimate of A<f>, and this best estimate may be corrected 
for slight variations in A<p at greater distances from the 
ground track. Such improvements would add to the 
complexity of the image processor, but would not de- 
part from the present invention. Consequently, it is 
intended that the claims be interpreted to include such 
modifications, improvements and equivalents as within 
the broader scope of the present invention. 

What is claimed is: 

1. A method for determining the elevation of terrain 
features in synthetic aperture radar image data using 
two antennae spaced apart a fixed distance, B, on a 
platform which carries the radar over terrain at a 
known altitude, a, with respect to a fixed reference, 
where the distance B is measured in a horizontal plane 
while the platform is in level flight, comprising the steps 
of 

separately processing the return signals from said 
antennae to the pixel level and determining the 
slant range, R, and phase difference, A<f>, of pixels 
from the return signals received through said two 
antennae from substantially the same point of said 
terrain, 

from said phase difference, A4>„ and said spacing, B, 
determining for each pixel the slant angle, 0, be- 
tween the slant range and a vertical line from the 
radar platform to the terrain, 

determining the vertical height, h, of said platform 
above each pixel from the slant range R as a func- 
tion of the cosine of the slant angle c f> for each pixel, 
and 

for each pixel, determining elevation, a', above said 
fixed reference as the difference between said 
known altitude a of the radar platform and the 
determined height h of said platform with respect 
to a corresponding point on said terrain for each 
pixel. 

2. A method as defined in claim 1 including the pro- 
cess of displaying said pixels with an indication of eleva- 
tion. 

3. A method as defined in claim 2 including the step 
of determining the horizontal distance, d, of each point 
on said terrain for each pixel from said vertical line and 
displaying each pixel at its distance from said vertical 
line. 

4. A method as defined in claim 1 wherein said indica- 
tion of elevation is displayed by quantizing elevation of 
pixels at predetermined constant intervals, and for pix- 
els of equal quantized elevation, displaying intensity in a 
predetermined color distinct from all other colors used 
for other quantized elevations over at least a distinct 
band of quantized elevations. 

5. A method as defined in claim 4 including the step 
of determining the horizontal distance, d, of each point 
on said terrain for each pixel from said vertical line and 
displaying each pixel at its distance from said vertical 
line. 
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6. A method for determining the elevation of terrain 
features in synthetic aperture radar image data using 
two antennae spaced apart a fixed distance, B, on a 
platform which carries the radar over terrain, where the 
distance B is measured in a horizontal plane while the 5 
platform is in level flight, comprising the step of 
separately processing the return signals from said 
antennae in complex vector form to the pixel level, 
and recording for each pixel from each antennae \q 
the altitude, a, of the radar platform above a fixed 
reference together with the roll angle, a, of the 
platform, 

determining the slant range, R, of pixels from the 
return signals received through said two antennae 15 
from substantially the same point of said terrain, 
from said pixel data in complex vector form, deter- 
mining the phase difference, A<j>, between pixels 
received from substantially the same point through 2Q 
said separate antennae from the relationship ex- 
pressed by the equation A<f>— (2irB/\) sin (0— a) 
where X is the wavelength of said radar, 
determining the vertical height, h, of said platform 
above each pixel from the slant range as a function 25 
of the cosine of the slant angle R for each pixel, and 
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for each pixel, determining elevation, a', with respect 
to said fixed reference as the difference between 
.said known altitude a of the radar platform and the 
determined height h of said platform with respect 
to a corresponding point on said terrain for each 
pixel. 

7. A method as defined in claim 6 including the pro- 
cess of displaying said pixels with an indication of eleva- 
tion. 

8. A method as defined in claim 6 including the step 
of determining the horizontal distance, d, of each point 
on said terrain for each pixel from said vertical line and 
displaying each pixel at its distance from said vertical 
line. 

9. A method as defined in claim 7 wherein said indica- 
tion of elevation is displayed by quantizing elevation of 
pixels at predetermined constant intervals, and for pix- 
els of equal quantized elevation, displaying intensity in a 
predetermined color distinct from all other colors used 
for other quantized elevations over at least a distinct 
band of quantized elevations. 

10. A method as defined in claim 9 including the step 
of determining the horizontal distance, d, of each point 
on said terrain for each pixel from said vertical line and 
displaying each pixel at its distance from said interval. 
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